Relaxation processes and entropic traps in the Backgammon model 
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We examine the density-density correlation function in a model recently proposed to study the 
effect of entropy barriers in glassy dynamics. We find that the relaxation proceeds in two steps 
with a fast beta process followed by alpha relaxation. The results are physically interpreted in the 
context of an adiabatic approximation which allows to separate the two processes, and to define 
an effective temperature in the off-equilibrium dynamics of the model. We investigate the behavior 
of the response function associated to the density, and find violations of the fluctuation dissipation 
theorem. 
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The relaxation in supercooled liquids near the glass 
transition has a characteristic two step form. Experi- 
ments on very different materials reveal the existence of 
a first fast relaxation process, called beta relaxation, fol- 
lowed by a much slower one, called alpha Jj]]. One of the 
most striking successes of the Mode Coupling Theory ||] 
is its ability to capture this phenomenon and to give a 
correct prediction for the relation among the exponents 
characterizing the two relaxations. However, we believe 
that a comprehension of the basic mechanisms underly- 
ing the relaxation in glasses is missing. Experiments in 
glasses have been recently interpreted in terms of traps 
models |3|,|]]. In these models, the system evolves among 
traps -or metastable states- which have a life time that 
grows decreasing the temperature, and finally diverges 
at the glass transition. The two step relaxation follows 
naturally from the hypothesis that equilibration inside 
a trap occurs much faster than "jumps" among differ- 
ent traps. How a trap may be defined and described for 
real systems or microscopic models is an interesting open 
problem. If the traps have to be interpreted as the result 
of energy barriers in a rough energy landscape, one finds 
the difficulty that the relaxation should appear as a ran- 
dom process even on a large scale. A "jump" among two 
different traps should imply a discontinuity in various 
quantities as the energy or the correlation function. This 
problem was already noted in [0], where it was proposed, 
as a way out, that real systems could be composed by a 
large number of quasi-independent sub-systems leading 
to the observed self averaging properties for the different 
quantities. In this direction, it can be instructive to in- 
vestigate a different mechanism for slow relaxation, and 
in particular the role of entropy barriers. 

In this letter we study the nature of density fluctu- 



ations in the so called "Backgammon model" , a micro- 
scopic model that has proved useful in studying some 
mechanisms underlying the glassy relaxation, and in par- 
ticular the role of entropy barriers. The model § is a 
Boltzmann gas with N particles in an N site space with 
Hamiltonian given by the total number of empty sites, 



JV 



H 



(i) 



where r — 1,...,N denotes the sites of the space and n r 
the occupation number of the site r. The system evolves 
following a single particle Metropolis dynamics: at each 
sweep a particle to move and an arrival site are chosen at 
random. The particle is moved with probability 1 if the 
energy does not increase, and with probability exp(— (3) 
if the move costs one unity of energy. The relaxation of 
the energy has been studied in detail in ]5| -j^], we briefly 
reassume here the main results. At zero temperature 
the dynamics is slower and slower as the times goes by: 
the average density of particles in the occupied states 
increases, and as a consequence the dynamics slows down. 
This observation has allowed for the identification of fast 
and slow degrees of freedom: the relaxation within the 
occupied sites at a given time proceeds much faster then 
the variation of the energy and the diffusion of "towers 
of particles" . This allows a self-consistent treatment of 
the dynamics which is in very good agreement with the 
exact solution. At large time the occupation probability 
P(n, t) mimics the equilibrium one 
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with effective time-dependent temperature T(t) = l/f3{t) 
larger than the effective one, and effective fugacity ((t) 
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related to (3(t) by the condition of constant density 
(n r (t)) = 1, i.e. + e«*) - 1 = CWe c(t) . At zero 

temperature, the decay of the energy follows the law: 
E(t) - -1 + 0(l/log(i)) while at small but finite tem- 
peratures this behavior is cutted-off for times of the or- 
der of the relaxation time r ~ exp{(3) / 1 when the re- 
laxation becomes exponential. Probing the system on 
finite time scales, e.g. mimicking heating cooling ex- 
periments, one finds characteristic glassy behavior, as 
hysteresis loops for the energy ||. Additional informa- 
tion on the off-equilibrium dynamics is gained studying 
the energy-energy autocorrelation function Cs(t, s) = 
<K^^>^m [|. This quantity shows aging 
at zero temperature with a scaling behavior Ce(£, s) ~ 
(s/t)^ (log(s)/ log(t)) ||. Again, at finite temperature 
this scaling is observed up to times of the order of the 
relaxation time. 

Here we concentrate our attention to the relaxation of 
the density-density correlation function, a quantity that 
is measured in experiments on real systems jjj . This is a 
better quantity to study fast processes in the system since 
on short time scales local densities varies while energy 
stays essentially constant. 

The fluctuations of the local density, around its average 
(n r (t)) = 1, are studied introducing the density-density 
correlation function 
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C«(M) = (n r (t)n r (s)) = 
n m P^ (n, t\m, s) (m, s) 



and its associated response 
i? w (i,s) = 



S(n r (t)) 
Sh r (s) 
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having denoted as P^ (n, t) the probability that state r is 
occupied by n particles at time t, P^ (n, t\m, s) the same 
probability conditioned to have m particles at time s, and 
h r is an infinitesimal inhomogcncous "pressure field" cou- 
pled linearly with the local density in the Hamiltonian 
pOf . We will present data for the normalized correlation 
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Choosing a site independent 



initial distribution, the correlation function remains site 
independent, while the response at finite time depends 
on the distribution on the h r i on the different sites, in 
the form of a dependence on = Prob(h r i > h r ) and 
v (r) _ p ro 5(/ lj ,, < /j r ) p Despite that, in the following 
we will drop the index r from the various quantities. At 
equilibrium correlation and response are time translation 
invariant and related by the fluctuation-dissipation theo- 
rem relation TR(t — s) = dC(t — s)/ds. In off-equilibrium 
conditions it has been often proved useful to characterize 
the violation of equilibrium by the "fluctuation dissipa- 
tion ratio" l|ll|,|ir 
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dC(t,s) 
ds 
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The evolution of the previously defined functions can 
be studied starting from the hierarchy of equations for 
P(n,t) and P(n,t\m, s) and related quantities following 
a procedure similar to the one used in . One can then 
derive closed integral equations in terms of few functions 
that can be integrated numerically or analytically in 
the long time limit [||. Details of this analysis, and 
in particular the full set of hierarchies for the density- 
density correlation (||) and the associated response func- 
tion (^) will be presented elsewhere. In this paper we 
integrate directly, truncating the hierarchy to some large 
order. In practice we have found that for low enough 
temperature and not too large times the truncation at 
n = 100 yields excellent results. 

Let us now discuss the form of the density correlation 
function in equilibrium. Starting from a random initial 
configuration, after times of the order of r Q ss e' 3 / 1 the 
system eventually reaches equilibrium. The correlation 
function is then time translation invariant and can be 
studied exactly in Laplace transform. It turns out that 
for temperatures small enough the dominant contribution 
to the equilibrium correlation function for all times is 
given by the superposition of two Debye processes 
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(t) = C 1 + ° e-^ 
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with r a = 2((C - l)e c + 1)/C 2 « 2e /3 //3 2 » C The fast 
and slow processes at equilibrium are related respectively 
to the relaxation within occupied states and to diffusion. 
We notice that the value of the plateau at equilibrium 
is (C — 1) / C; that corresponds to complete decorrelation 
within occupied states, but no diffusion. The equilibrium 
curve for T = 0.05 is the dashed line in fig. 1. In the 
time window shown in figure 1 the correlation function 
stays essentially constant equal to the plateau value (£ — 
1)/C ~ 0.94 decaying to zero afterwards. Here the alpha 
relaxation time is r Q ~ 3 x 10 6 . 

Regarding the off-equilibrium correlation function we 
show in figure 1 the data for the C norm (t, s) for T = 0.05 
as a function of t — s for various values of s. Over the 
time window we explore the system is far from being 
thermalized. The integration at T — on the same time 
window leads to almost identical results. The a and (3 
processes are well separated for large enough s. In this 
case, the shape of the off-equilibrium relaxation curve 
can be understood qualitatively within the framework 
of the adiabatic approximation. The decorrelation time 
for the density among filled states is much smaller that 
the time needed to diffuse and/or change sensibly the 
energy; there must be then a time scale such that we can 
approximate 
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P(n,t\m,s) w S mfi S nfi + (1 - 5 m0 )(l - S n ,o) 



and consequently 
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l--Po(s) 
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For short times (t — s « s), ((s) = — log(Pi(s)) does 
not vary too much; it is then natural to describe the beta 
relaxation by the form (^|) just substituting £ by £( s )- 
In this way we have checked we correctly describe even 
the beginning of the alpha relaxation. However the best 
combined description of the slow a relaxation together 
with the [3 relaxation is obtained by functions of the aging 
form: 
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C(s) -l + e-(*-">/<W (1 + &00) 
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a form inspired by the one of the energy-energy correla- 
tion function found in At finite temperature, b(s) is 
a crossover function to the form (^|) . 
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FIG. 1. Density-density correlation function C n orm(t, s) as 
a function of t — s for different values of s = 1, 10, 10 2 , 10 3 , 10 4 
(from left to right) at T = 0.05 (continuous lines). The exis- 
tence of a plateau separates the beta and alpha regimes. The 
dashed line corresponds to the equilibrium C^o rm (t — s). The 
curves with s = 10 3 ,10 4 are excellently fitted by the form 
(0) with, respectively, b(s) = 2.20, 3.85. For the same times, 
C(s) = - log P(l, s) = 8.93, 11.51. 

In figure 2 we see that the two relaxation processes 
are also manifest in the response function. In order to 
characterize better the two processes in off-equilibrium 
conditions we study the relation among the response and 
the correlation during the dynamics. The first quantity of 
interest is the value of the fluctuation-dissipation ratio at 
equal times. We see that x(t, t) reaches values close to one 



much before total equilibrium sets in. This is a further 
indication that the system is in local equilibrium and 
yields in a natural way the notion of a trap in this system. 
Note that the traps in the backgammon model are purely 
cntropic. This means that the system escapes from the 
trap even at zero temperature when thermal excitations 
are absent. In other words, the dynamics is slowed down 
by entropic traps which can be considered as metastable 
states with a finite (energy dependent) lifetime. 
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FIG. 2. The response function for the values of 
s = 1, 10, 10 2 , 10 3 , 10 4 (from top to bottom). 
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FIG. 3. Fluctuation-dissipation ratio x(t,s) for different 
values of s = 10, 10 2 , 10 3 , 10 4 (from bottom to top) as a func- 
tion of 1 — C n0 rm(t, s) . The two plateaux correspond to the 
alpha and beta relaxations. 

Most interestingly we find that for large times (large 
values of s) x(t, s) remains nearly constant in each of the 
two processes. This is shown in fig. 3 where we plot the 
x(t, s) as a function of 1 — C norm (t, s) for different values 
of s. We find that the x(t, s) is close to 1 in the beta 
relaxational process and jumps to a s-dependent smaller 



3 



value ^a(s) at a value Cnormi,^: ■->) — Q_ea remaining con- 
stant in the alpha process. 

It is rather natural, with the adiabatic approximation 
in mind, to look for a possible interpretation of x as the 
ratio between the effective temperature and the actual 
one. Unfortunately we found negative evidence for such 
interpretation of the data. For low enough temperature 
and times smaller than t q the actual value of x(t, s) is 
nearly temperature independent. The step behavior of 
the x(t, s) is quite reminiscent of the dynamics in spin 
glasses with one step of replica symmetry breaking [ p"T| . 
This in turn represents the suitable off-equilibrium gener- 
alization of Mode Coupling Theory E3] in case of Whit- 
ney fold glass singularity, i.e. precisely the case where 
there are well separated a and f3 relaxations. We do not 
know at this stage if this corresponds to a deep analogy 
or is a mere coincidence. We stress that at equilibrium, 
the ending of the beta relaxation, as well as the onset 
of the a relaxation are described by power laws in Mode 
Coupling Theory. As it is seen in (|^), in our case the 
equilibrium relaxation is the superposition of two expo- 
nential. 

Summarizing, we have shown the existence of two step 
relaxational processes (alpha and beta relaxation) in the 
Backgammon model, which is a simple model where the 
slow dynamics is consequence of pure entropic barriers. 
We have closed the dynamical equations for the density- 
density correlation function and the response of the sys- 
tem to a staggered field coupled to the density. Our 
results for those quantities and the fluctuation dissipa- 
tion ratio allow for a clear identification of the fast (beta) 
and slow (alpha) relaxation processes in this system. The 
physical interpretation of this processes leads naturally 
to the concept of entropic trap. In case of entropic traps, 
the entropy barrier associated to the trap depends it- 
self on the number of available configurations within the 
trap. Because the height of the entropic barriers varies it- 
self in a continuous way, the dynamics proceeds slow but 
without macroscopic jumps in the energy. The results 
presented here for the Backgammon model are expected 
to apply for slowly relaxing systems where the dynamics 
is mainly driven by entropy barriers (like for instance, 
Bose-Einstein condensation) . A detailed account of our 
work will be given elsewhere. 
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